Breakup of the Fermi surface near the Mott transition in low-dimensional systems 
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We investigate the Mott transition in weakly- coupled one-dimensional (Id) fermionic chains. Us- 
ing a generalization of Dynamic Mean Field Theory, we show that the Mott gap is suppressed at 
some critical hopping t^ 2 . The transition from the Id insulator to a 2d metal proceeds through an 
intermediate phase where the Fermi surface is broken into electron and hole pockets. The quasipar- 
ticle spectral weight is strongly anisotropic along the Fermi surface, both in the intermediate and 
metallic phases. We argue that such pockets would look like 'arcs' in photoemission experiments. 



The formation of a pseudogap phase near the Mott 
transition in strongly correlated electron system is a long- 
standing question in condensed matter physics. Since the 
observation of a pseudogap state in high-T c copper-oxide 
superconductors by a variety of experimental techniques 
^J, the nature and origin of this phase lies at the core 
of a heated debate. Several related issues, such as the 
possible breakup of the Fermi surface (FS) into arcs or 
pockets 0, Q and the concomitant appearance of 'hot 
spots' in the quasiparticle spectrum, have also attracted 
much attention in recent years. The theoretical descrip- 
tion of the Mott transition in two spatial dimensions has 
proven to be difficult |4|, however, and a consensus on 
the nature and properties of the pseudogap has not yet 
emerged. The central question concerns the way the FS 
is destroyed when entering the Mott state, either at zero 
temperature as a function of some control parameter, 
or with decreasing temperature from a high-T metallic 
phase. 

For solving these questions one has to resort to ap- 
proximate or numerical approaches. Among those meth- 
ods, the Dynamic Mean Field Theory (DMFT) HQ has 
proven to be very fruitful to tackle the question of Mott 
transitions and strong correlations in high dimensional 
systems. However, since in the standard implementa- 
tion of this method one reduces the problem to a single 
site, the momentum dependence of the self energy is lost, 
which makes it unable to describe the FS anisotropics, 
and answer the above questions. To overcome this lim- 
itation, various cluster-DMFT schemes have been used 
to approach the 2d problem 0, H, E3, [HJ . One difficulty 
with these approaches which, for small clusters, can af- 
fect physical predictions, lies in a certain degree of ar- 
bitrariness in converting the cluster results into physical 
quantities on the lattice. 

In addition to high-T c superconductors, these ques- 
tions related to the approach of the Mott transition are 
directly relevant to other systems in which interactions 
are deemed to be important such as the organic supercon- 
ductors, made of weakly-coupled one-dimensional chains 



[lH . Such systems exhibit a deconfinement transition 
between a one-dimensional Mott insulator and a three 
dimensional metal. For such systems of coupled chains, 
FS pockets were found when the inter-chain coupling 
is treated within the RPA approximation 12]. In that 
model the pockets form when the inter-chain bandwidth 
exceeds the Id Mott gap of the chains. However the RPA 
approximation neglects the feedback of the inter-chain 
hopping on the Id Mott gap, and more generally on the 
self-energy. Depending on the evolution of the Mott gap 
with inter-chain coupling, the FS pockets may or may 
not form. To go beyond this approximation, we have 
proposed a generalization of the one-site DMFT to one- 
dimensional chains embedded in a self-consistent bath 
0. This ch-DMFT scheme was successfully applied to 
study the deconfinement transition in quasi- Id fermionic 
lattices, and is the ideal tool to investi gate the Mott tran- 
sition in this type of systems 0, fT^L fT^| . This method 
has the advantage of having no arbitrariness in the def- 
inition of the cluster, since here the cluster corresponds 
to a physically well-defined entity (a single chain). 

In the present paper, we use this method to study 
a two-dimensional lattice of spinless fermions made of 
weakly coupled one-dimensional chains. Contrarily to 
the Id Hubbard model, the Id spinless model has a Mott 
transition for a finite value of the interaction. Further- 
more, the effects of the interaction are stronger in the 
spinless case than in the spin--| case, as illustrated by 
the lower value of the Luttinger coefficient which can 
be reached in the former. This makes the spinless model 
very well suited to tackle these issues. The spinless model 
can also be viewed as a caricature of a spin-| model at 
quarter filling with a very strong local repulsion. Using 
the ch-DMFT method we show that FS pockets exist in 
this model close to the metal-insulator transition. How- 
ever, due to the very strong anisotropy of the spectral 
function and of the spectral weight along the Fermi sur- 
face, these pockets would appear as Fermi-surface arcs in 
experiments. We note that Stanescu and Kotliar [13 and 
Kuchinskii et. al. 18] have discussed similar effects using 
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extensions of DMFT including momentum dependence. 

Our model has two parameters, t± and V, which con- 
trol the inter-chain hopping and the interaction 7Yi n t = 
V n r n r+ i within the chains, respectively, relative to 
the in-chain hopping t. In the Id limit (t± = 0) this 
model has a phase transition between a Luttinger liquid 
at V < 2 1 and a Mott insulator at V > 2t. Here we focus 
on the region V > 2 1, and we investigate the destruction 
of the Mott insulator with increasing t±. Within ch- 
DMFT the 2d problem is mapped onto an effective Id 
problem described by the action 
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The inverse propagator £f -1 in Eq. plays the role 
of a long-range, time-dependent hopping amplitude. It 
must be determined from the requirement that the in- 
chain Green's function calculated from S e ff coincides 
with the k \ -summed Green's function of the original 2d 
model |3,LL4(- This requirement implies 
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with k the in-chain momentum, ^ 
bare dispersion, and R(z) = sign[Re(z)] y/l/z 2 + (2t± 



resulting from the integration over the transverse energy 
e± = — 2t± cos(k±). Eqs. JJJ and (|2j) can be readily 
derived from the assumption that the lattice self-energy 
does not depend on transverse momentum k±, and is 
given by the effective in-chain self-energy: 



E(fc,u;) = < ^ ~ 1 (k,co) - < ^~ 1 (k,cj). 



(3) 



E(fc,a;) will be our main concern here. In order to eval- 
uate the self-energy we compute the space-time Green's 
function Sf(r, r) = — (c r (r)cj(0)) of the action Eq. (0) by 
quantum Monte- Carlo on a discrete imaginary-time mesh 
T£ = £(3/L, using the Hirsch-Fye algorithm From the 
Fourier transform £f (fc, ico n ) we construct a new propaga- 
tor % according to Eq. (|2j) , which is fed back into Eq. JU 
until self-consistency is achieved. 

We consider a half-filled 32-sites chain closed with an- 
tiperiodic boundary conditions. These boundary condi- 
tions were adopted for two reasons: (i) in short Id chains 
we observed that the convergence of the Monte-Carlo to- 
ward exact-diagonalization results is much faster with 
antiperiodic than with periodic boundary conditions; (ii) 
for a given system size, the antiperiodic boundary con- 
ditions improve the resolution near k = 7r/2, which is 
an advantage for the investigation of the FS properties. 
Indeed, the shape of the FS is controlled by the real part 
of the self-energy through the equation 
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FIG. 1: Left panels: Self-energy on the imaginary- frequency 
axis for k between and tt/2 at V = 4, t± = 0.5, and T = 
0.04. The dots show the numerical data and the lines are 
guides to the eye. The number of imaginary-time slices was 
L — 60. Right panels: Fit of the numerical data to a trial 
self-energy (see text). All energies are in units of t. 



which must be solved at zero temperature near k = tt/2. 
At finite temperature the FS looses its identity, although 
sharp signatures may subsist in the zero-energy spectral 
function A(k, k±, uj = 0), which is accessible through 
photoemission experiments. We will first discuss the FS 
topology and properties implied by Eq. (@J and our ch- 
DMFT results for E(fc,a;), before addressing some issues 
related to the experimental measurement of the FS. 

For evaluating Eq. one needs a procedure to con- 
tinue the ch-DMFT self-energy from the lowest Matsub- 
ara frequency zujq = iirT down to uo = i0 + along the 
imaginary- frequency axis. This is a delicate endeavor, 
which in general requires some assumption about the 
analytical behavior of E(fc,o;) near uj = 0. We per- 
formed the continuation by fitting the self-energy to 
an analytical function. The prominent feature in the 
Mott phase is the spectral gap which can be crudely 
represented by a self-energy of the form Ea(&,c<;) = 
(A - + (A/2) 2 /(cj + A&), where A is the gap width 
and A accounts for the dispersion renormalization due 
to exchange. A similar Ansatz was recently proposed to 
describe high-T c cuprate superconductors [2(| . This sim- 
ple form is not sufficient to reproduce our QMC results, 
however, even in the pure Mott phase at t± = 0. We 
obtained a much better agreement with our data by tak- 
ing into account residual interactions. Specifically, the 
self-energy to which we fit the QMC results takes the 
form E(fc,o;) = Ea(&,c<;) + Ei n t(fc, cj), where Ei nt con- 
tains all diagrams — evaluated using the gapped electron 
propagator Go(k,uj) = [uj — — Ea(&,^)] -1 and an ef- 
fective interaction strength V* — up to second order in 
perturbation theory [2]] . A comparison of the QMC and 
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model self-energies is displayed in Fig. [I] One can see 
that the model has enough freedom to fit the QMC re- 
sults in great detail, especially in the low-frequency re- 
gion we are mostly interested in. Furthermore, it turns 
out that the model fits the QMC data in the whole range 
of temperature and t± values which we have investigated. 
We can therefore use this fit to track the closing of the 
Mott gap and the formation of the Fermi surface as t± 
is increased. A few additional illustrations of the fit per- 
formance can be seen in Fig. 

In Fig. [21 we display our results for — Re E(fc, icoo) at 
V = 2.5 1 and different values of t±, together with the fit 
results evaluated at icoo and i0 + . The main trend with 
increasing t± can be seen on the raw numerical data. 
For small t± the self-energy has a tendency to diverge 
near fcp = ^/2 (Fig. Et). This behavior is most clearly 
seen at large V and low T (see Fig. Q). The singularity 
of ReE(fe, 0) results from the presence of a gap in the 
zero-temperature spectral function at fcp E3 •> an d is well 
captured by the model self-energy evaluated at T = 
(blue line). At finite frequency and/or temperature the 
singularity is regularized as shown by the red line. With 
increasing t± the drop of the QMC self-energy across kp 
diminishes (Fig.Eb). Correspondingly the fitted spectral 
gap A decreases and eventually vanishes at tj 2 ~ 0.5 1, 
together with the disappearance of the singularity in the 
self-energy (Fig. Efc)- 

Fig. [21 provides the graphical solution of Eq. and 
illustrates the formation of the Fermi-surface pockets in 
this model. When t± is small, the transverse dispersion 
is not sufficient to overcome the gap in the self-energy, 
and the system remains insulating. At high on the 
contrary, there is no gap in the self-energy and Eq. (0} 
has a solution for all fcj_, leading to a continuous Fermi 
surface. The latter has practically the same shape as 
the non-interacting FS, but is strongly renormalized to 
an effective inter-chain hopping tj_ ~ 0.41 1±. In the 
intermediate regime t± < t± < t^ 2 , there is a finite range 
of k± values where Eq. (gj) admits two solutions, leading 
to the breakup of the Fermi-surface into pockets (Fig. 03 
and b'). 

According to Luttinger's theorem, the area of the Bril- 
louin zone where Re£f (fc, k \ , 0) > equals the electron 
density and is thus conserved [lj, [23| . In the Mott phase 
Re£f (fc, fc_L, 0) changes sign at k = ±tt/2 due to the di- 
vergence of ReE(fc,0) (Fig. Ek), and is positive in the 
domain \k\ < 7r/2, leading to a density n = 1/2 as ex- 
pected. Because the singularity of ReE(±7r/2,0) sub- 
sists as long as A > 0, the line of zeros of Re£f(fc,0) 
at k = ±7r/2 is still present when FS pockets develop, 
as indicated in Fig. On the other hand, owing to 

particle-hole symmetry — which reduces at uo = to an 
inversion about the point (fc, k±) = (tt/2, tt/2) — the elec- 
tron and hole pockets have identical volumes, so that 
Luttinger theorem is trivially obeyed in our results. 

The mechanisms of FS pockets formation in the present 
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FIG. 2: (a-c) Real part of the ch-DMFT self-energy as a 
function of k at the lowest Matsubara frequency for V = 2.5, 
T — 0.1, and increasing t± (red points). The red and blue 
lines show the fit of the data to the function Ea(&,^o) + 
Eint(fc, iujq), and its continuation to zero energy, respectively. 
For the fit the model self-energy was evaluated on a 32 k- 
point mesh, and then reevaluated using the fitted parameters 
on a dense mesh at uj = iuo and T — 0.1 (red lines) and 
at co = i0 + and T = (blue lines). The shaded areas show 
the domain covered by the free dispersion — 2t± cos(k±) in 
Eq. @. The Fermi surface corresponding to (b) and (c) is 
shown in (b') and (c'), respectively. The dotted line indicates 
the non-interacting Fermi surface in each case, and the + and 
— indicate the sign of ReSf k±, 0) within the Brillouin zone. 



study and in the RPA approach of Ref. are similar, 
although there is one important difference. In the RPA 
approximation the spectral gap keeps its Id value at all 
tj_: pockets form when t± > t c ^ ~ Aid, and they never 
merge into a connected Fermi surface as t± continues 
increasing. Within ch-DMFT, in contrast, the closing of 
the Mott gap with increasing t± is correctly captured; 
as a result the pockets form at lower t± values — thus 
they are very thin — and eventually they disappear at tj 2 
where A = 0. 

The quasiparticle properties are strongly anisotropic 
along the Fermi-surface pockets. It is already clear from 
Fig. [213 that the spectral weight is much smaller on the 
vertical parts of the pockets closest to k = 7r/2, due to 
the diverging self-energy in this region; as a result the 
pockets would most likely look like 'arcs' in photoemis- 
sion experiments (see below). Fig.|3Ji shows the evolution 
of the quasiparticle residue along the Fermi surface. The 
residue was evaluated as Z = [1 — dRe E(fe, ou) / duo^^) -1 
using the model self-energy and an interpolation of the 
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FIG. 3: (a) Evolution of the quasiparticle residue Z along 
the Fermi surface for different inter-chain couplings t±. When 
FS pockets are present there are two values of Z for each k± , 
the lowest value corresponding to the region of the pocket 
closest to k = 7r/2. (b) Parameters of the model self-energy 
determined from fits to the ch-DMFT numerical data. 
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FIG. 4: Comparison of the zero-energy spectral function 
A(k,0) at T = with the expected ARPES intensity I(k) 
calculated assuming a /c-space resolution of 0.047T, an en- 
ergy resolution 0.004 1, and an energy integration window 
5E — 0.01 1. The self-energy parameters are taken from 
Fig. Eb for t± = 0.42 (a) and t± = 0.46 (b). 



parameters fitted to the ch-DMFT data in the whole 
range of t± values (Fig. Eb). In the intermediate phase 
the residue on the vertical parts of the pockets decreases 
from Z ~ 0.15 to Z ~ with increasing t±. On the 
'cold' side of the pockets the behavior is inverted, and the 
residue increases from Z ~ 0.35 to Z ~ 0.5. Strikingly, 
a hot spot around k = k± = tt/2 subsists at t± > tj 2 . 
Here again the evolution of Z with t± is different at the 
cold and hot spots: while Z steadily approaches one in 
the cold region, it remains close to Z ~ 0.4 at the hot 
spots. 

It's worth stressing the role of residual interactions in 
the results of Figs. and OH At the qualitative level, 
the self-energy Sa(^,w) (with A = 1) is sufficient to un- 
derstand the formation of FS pockets with anisotropic 
residues. Indeed, using E = Ea in Eq. (gj) one finds 
that pockets form for < A < 2tj_, and that the residue 
varies on such pockets. However the pockets obtained in 
this way are considerably wider than in Fig. [2J and the 
residue in the cold regions is Z ~ 1 when A approaches 
zero instead of Z ~ 0.5 as in Fig. 03 Thus the residual in- 
teractions are important for the quantitative understand- 
ing of the FS properties. Meanwhile, the fact that our 
model self-energy fits the ch-DMFT data at t± > t*f with 
V* ~ V (Fig. Eb) indicates that second-order perturba- 
tion theory is a good approximation in this region, as 
expected in a Fermi liquid. 

We now turn to the question of the experimental ob- 
servation of FS pockets. There are several limitations 
which could make the observation of such pockets by 
ARPES challenging, such as the finite energy and mo- 
mentum resolutions, the finite temperature at which ex- 
periments are performed, as well as the necessity to inte- 
grate the ARPES intensity on some energy window in or- 
der to improve the signal to noise ratio. Ideally, ARPES 
would measure the occupied spectrum A(k,uj)f(uj). In 
practice, however, due to the above limitations, the mea- 
sured intensity at the Fermi energy would be I(k) oc 



J^° §E du; J dedq A(q,e)f(e)g(k — q,u — e) where SE de- 
fines the energy integration window and g is some func- 
tion describing the instrument resolution. We have cal- 
culated I(k,k±) using a Gaussian for g. The comparison 
depicted in Fig. 0] of the T = Fermi surface with the 
expected ARPES intensity clearly shows that the closing 
segments of the pockets near k = n/2 would very likely 
be hidden in the ARPES signal. The broad aspect of 
J(fe) as compared to A(k,0) is not a consequence of fi- 
nite temperature, but of (i) the finite /c-space resolution 
combined with the fact that the pockets are very thin 
and (ii) the large difference in quasiparticle weight on 
the two sides of the pockets, which is obvious in A(k, 0) 
and consistent with the residues shown in Fig. [3 Simi- 
lar FS anisotropics were recently found in cluster-DMFT 
studies of the 2d Hubbard model 0, suggesting that 
such effects are generic to systems close to a Mott transi- 
tion, and could possibly explain the ARPES observation 
of FS arcs in high-T c cuprate superconductors. 
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